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/  A  new  algorithm  for  computing  scattering  kernels  from  Arctic  ice 
keels  is  described.  The-  basic  equations  are  adapted  from  Burke  and 
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Section  1 
INTRODUCTION 


The  Arctic  ice-covered  environment  presents  a  sort  of  surface 
scattering  that  is  quite  different  from  that  in  the  open  ocean.  Under-ice 
acoustic  propagation  undergoes  incoherent  scattering  in  the  presence  of 
ice  keels.  Greene  and  Bowen  (1983)  reviewed  the  special  environmental 
conditions  that  differentiate  acoustic  propagation  in  the  Arctic  from  open 
ocean  regions.  The  under-ice  roughness  associated  with  ice  keels  has  a 
very  important  effect  on  sound  propagation.  Greene  (1984)  developed  a 
statistical  model  of  under-ice  roughness,  based  on  a  distribution  of 
randomly  oriented  cylindrical  keels.  Using  this  model,  Greene  developed 
an  interim  scattering  model  (SISM/ICE),  which  is  a  hjtorid  of  theories 
based  on  continuous  and  discrete  roughness  models.  This  type  of  charac¬ 
terization  of  scattering  is  important  for  the  Navy's  propagation  models. 
Rubenstein  et  al.  (1986)  described  an  implementation  of  SISM/ICE  for  the 
ASTRAL  and  PE  models.  The  concepts  are  also  relevant  to  other  models, 
including  FACT,  FFP,  MPP,  and  RAYMODE. 

The  SISM/ICE  model  is  probably  adequate  for  low  frequencies, 
below  about  200  Hz.  However,  at  higher  frequencies,  SISM/ICE  ignores  the 
scattered,  incoherent  component  of  energy  that  is  directed  within  the 
sound  channel  (Greene  and  Rubenstein,  1986).  Therefore,  a  more  realistic 
scattering  kernel  is  needed  for  higher  frequency  applications.  In  this 
report,  an  algorithm  is  described  for  computing  realistic  scattering 
kernels  for  a  field  of  ice  keels. 

The  under-ice  surface  is  considered  to  be  flat,  with  cylindrical 
bosses  of  elliptical  cross  section.  The  basic  assumptions  are  that  the 
ice  keels 


(1)  have  a  constant  size,  and  a  half-width-to-depth  ratio  of 
1.6,  as  suggested  by  Diachok  (1976), 
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(2)  have  random  orientation, 


(3)  have  random  spacing  with  uniform  distribution  along  a 
track,  and  that 

(4)  pressure  release  boundary  conditions  are  appropriate. 

Section  2  of  this  report  describes  the  model  equations  in  some 
detail.  The  basic  equations  are  from  Burke  and  Twersky  (1966).  A  new 
method  of  calculation  is  described  in  Section  3.  The  method  involves  a 
direct  calculation  of  power  series  of  Mathieu  functions,  and  is  more 
accurate  than  the  approximate  methods  of  Burke  and  Twersky.  Section  4 
presents  results  of  applying  the  model  to  determine  reflection  coeffi¬ 
cients,  scattering  cross  sections,  and  scattering  losses  for  particular 
fields  of  ice  keels. 
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Section  2 
MODKL  EQUATIONS 


2.1  MATHIBU'S  BQOATXOMS  AMD  MATBZBD  PUECTIOHS 

The  elliptic  coordinate  system  is  a  natural  one  to  use  for  the 
problem  of  scattering  off  an  elliptical  cylinder.  Following  the  notation 
from  Morse  and  Feshbach  (1  953),  the  coordinates  y  and  0  are  related  to 
Cartesian  coordinates  by 

x  *  1 /2  a  cosh  y  cos0, 

y  *  1/2  a  sinh  y  sin0,  (2.1) 

where  a  curve  y  *  constant  is  an  ellipse  and  a  curve  0  =  constant  is  a 
hyperbola,  each  with  foci  at  x  =  ±  1/2  a.  See  Figure  2.1  for  a  schematic 
drawing.  For  reference,  the  elliptic  coordinates  are  related  to  polar 
coordinates  by 

r  =  1/2  a  [cosh2  y  -  sin20]1/2  , 

$  *  tan-1  [tanh  y  tan©]  .  (2.2) 

In  Cartesian  coordinates,  the  Helmholtz  wave  equation  is 

written 


(¥-  ♦  ♦  k2)  *  *  0  ,  (2.3) 

3x2  3y2 


and  transforming  to  elliptic  coordinates,  the  equation  becomes 

lit  +  h2  ( cosh2y  “  cos20)  iji  =  0  ,  (2.4) 

3u2  392 


where 


k  *  2ir/X  is  the  wavenumber,  and  X  is  wavelength.  Bq.  (2.4)  separates  into 
two  ordinary  differential  equations 

j2u 

+  (b  -  h2cos20)H  -  0  ,  (2.6) 

d92 

j2u 

-  +  (b  -  h2cosh2u)M  *  0,  (2.7) 

du2 


with 


<|>  *  H( 0 )  M(u)  •  (2.8) 

Bq.  (2.6)  is  Mathieu’s  equation,  and  (2.7)  is  known  as  Hathieu's  modified 
equation. 


The  solutions  to  (2.6)  are  periodic  in  0.  The  solutions 

00 

Se  (h,0)  *  y  A  (h,m)  cos  (n0)  , 

n  „  n 

n-0 

00 

So  (h,0 )  *  V  B  (h,m)  sin  (n8)  ,  (2.9) 

m  .  n 

n=1 

are  even  and  odd  functions,  respectively,  about  6*0.  Those  functions  with 
even  order  m  have  period  n,  and  those  with  odd  order  m  have  period  2ir. 
Also,  the  even  order  functions  are  even  about  0  *  v/2  and  3ir/2,  and  the 
odd  order  functions  are  odd  there. 


The  Fourier  coefficients  in  (2.9)  are  normalized  as 


In  addition,  normalization  constants  are  defined  by 


He  =*  Je  +  i  Ne  , 
n  n  n 


Ho  *  Jo  +  i  No  . 
n  n  n 


SCATTERING  SOIOTICM 


(2.13) 


Ne  develop  the  solution  for  plane  waves  scattered  by  a  semi¬ 
elliptical  cylinder.  In  terms  of  Mathieu  functions,  a  plane  wave  is 
expressed  by 

ikr  cos(u-a)  ik  (xcosu  +  ysinu) 


»  Sj  (h,u)Sj  (h, 0) 

&  l  l  *  <*'»>' 


n*0  j»e,o 


Mjn(h) 


(2.14) 


where  the  summation  variable  j  indicates  a  summation  over  even  and  odd 
Mathieu  functions,  and  the  subscript  i  denotes  a  wave  that  is  incident  on 
a  cylindrical  surface.* 


The  scattered  wave  (see  Figure  2.2)  is  represented  in  terms 
of  outgoing  radial  Mathieu  functions, 

__  «  Sj ' (h,u)Sj  (h,0) 

l  .  I  wThT -  HV1”“)  '  (2',5) 

n=»0  :=e,o  n 

where  the  coefficients  Xj  are  determined  by  the  boundary  conditions.  We 

n 

use  a  free  surface  boundary  condition, 


h  +  *s  * 0 


(2.16) 


*  This  expression  for  <j^  is  the  same  as  that  quoted  by  Morse  and  Feshbach 
(1953),  but  is  different  from  that  quoted  by  Burke  and  Twersky  (1964) 
by  a  factor  /  tt/ 2 .  However,  a  second  discrepancy  in  Burke  and  Twersky 
(1964)  in  their  asymptotic  expression  for  radial  Mathieu  functions  Hem 
and  Hom  compensates  exactly,  so  that  the  later  results  fall  into  agree¬ 
ment. 


i v/ aV-V  j  X  A' t  .Vj\*  o  .*_•  ‘j>.  \*\V  *.»  njA -V .  -  .  • .  ■ 


Diagraa  of  elliptical  cross  section  of  cylinder  in  the  lower 
half-plane,  incoming  grating  angle  u  and  outgoing  angle 
The  angle  -u  corresponds  to  virtual  image  plane  wave. 


Therefore,  the  coefficients  become 


JJ„<h,u0) 

X3n  "  '  Hj  (h,u  )  ' 
n  o 

where  u  *  u  is  the  surface  of  the  elliptical  cylinder, 
o 


(2.17) 


To  obtain  the  solution  for  a  semi-elliptical  cylinder,  we  use 
Rayleigh's  image  technique.  we  superpose  two  incident  waves  from  direc¬ 
tions  u  and  -u,  and  compute  the  total  scattered  field.  Because  of  the 
symmetry  properties  of  periodic  Mathieu  functions, 

Se  (u)  »  Se  (-u)  , 
n  n 


So  (u)  =■  -So  (-u)  , 
n  n 


(2.18) 


the  total  scattered  field  becomes 


T  ■  ♦  (u)  -  iL  (-u) 
s  s 

_  •  So  (h,u)  So  (h,9) 

-  2  /8*  l  inXo  -n  -  " -  Ho  (h,u)  .  (2.19) 

n-0  n  %ih)  n 

We  can  derive  an  approximation  to  (2.19)  which  is  accurate  in 
the  far  field.  Using  the  asymptotic  formulas 

So  (h,0)  So  (h,a)  , 
n  n 

Ho  (h,y)  *  ✓  -r~—  eilcr  i"n  ,  (2.20) 

n  lKt 


we  write 


_ 2_  ikr 

f  =  /  iwkr  e 


f  (  $,u) , 


(2.21 ) 
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where 


co  So  (h,u)  So  (h,$) 

(  f(*,u)  -  4*  I  Xon  MO  (h)n 

I  n=0  n 

t 


(2.22) 


i  We  have  written  the  function  f(^,u)  in  a  form  comparable  to  the  f_  func¬ 

tion  used  by  Burke  and  Twersky  (1964,  1966),  with  u  and  ^  being  the  incom¬ 
ing  and  outgoing  grazing  angles,  respectively.  The  normalization  con¬ 
stants  Mo  (h)  are  defined  as  in  Morse  and  Feshbach  (1953)  and  Burke  and 
n 

Twersky  (1966),  but  are  different  from  those  used  by  Burke  and  Twersky 
(1964)  by  a  factor  2ir.  We  also  note  that  our  angles  here  are  grazing 
angles,  whereas  in  Burke  and  Twersky  (1964,  1966)  they  are  incident. 


We  call  the  case  treated  here  oblate,  because  the  ground  plane 
y»0  is  aligned  parallel  to  the  elliptical  cylinder's  major  axis.  Because 
of  the  symmetry  of  this  case,  the  even  Mathieu  functions  have  dropped  out 
of  (2.19)  and  (2.22).  The  prolate  case,  in  which  the  major  axis  is 
aligned  perpendicular  to  the  ground  plane,  is  also  treated  by  Burke  and 
Twersky  (1964,  1966).  Because  it  is  not  specifically  of  interest  to  the 
ice  keel  scattering  problem,  it  is  not  included  here.  The  expressions 
analogous  to  (2.19)  and  (2.22)  for  the  prolate  case  involve  summations 
over  odd-order  even  periodic  Mathieu  functions  and  even-order  odd  periodic 
Mathieu  functions. 

2.3  GRID  OF  RANDOMLY  DISTRIBUTED  SCATTBUERS 


Here  the  equations  for  a  grid  of  randomly  distributed  ice  keels 
(semi-elliptical  cylinders)  are  summarized,  from  Burke  and  Twersky  (1966). 
We  define  a  quantity 


Z 


n 

k  sinu 


f(u,u)  , 


(2.23) 


where  n  is  the  number  density  of  scatterers  per  unit  length.  Then  the 
coherent  power  reflection  coefficient  is  given  by 


v  *:«, 


i  +z  i  2 


The  incoherent  scattering  cross  section  per  unit  length  is  written 


(2.24) 


2n  .fU,u) |2 
0(*'u)  "  * 


Conservation  of  energy  may  be  expressed  by 


R  +  csc(u)  /  <j(  u)  d$  ■  1  , 


(2.25) 


(2.26) 


that  is,  coherent  reflected  energy,  plus  the  integral  of  incoherent 
scattered  energy,  equals  the  total  incident  energy. 


Section  3 

unoo  or  calculation 


Burke  and  Twersky  (1966)  used  approximate  methods  for  computing 
the  scattering  solutions,  their  approximations  are  valid  in  the  limits  of 
low  and  high  frequencies. 

Instead  of  using  these  approximate  methods,  we  calculated  the 
reflection  coefficient  R  and  scattering  cross  section  o  directly,  by 
explicitly  evaluating  the  power  series  (2.22). 

The  algorithm  developed  by  Clemm  (1969)  was  used  to  compute  the 
periodic  and  radial  Mathieu  functions,  this  algorithm  uses  trigonometric 
function  power  series  for  computing  periodic  Mathieu  functions,  and  Bessel 
function  power  series  for  computing  radial  Mathieu  functions.  In  comput¬ 
ing  these  power  series,  a  tolerance  level  of  10“13  is  set.  the  power 
series  is  truncated  when  the  magnitudes  of  two  successive  terms,  relative 
to  the  magnitude  of  the  largest  term,  is  less  than  or  equal  to  the  toler¬ 
ance  level.  Test  cases  were  run,  comparing  computed  values  of  Mathieu 
functions  with  values  tabulated  in  Abramowitz  and  Stegun  (1964).  Agree¬ 
ment  was  good  to  seven  decimal  places.  In  order  to  achieve  agreement  with 
Morse  and  Feshbach  (1953),  it  is  necessary  to  choose  the  Ince  normaliza¬ 
tion  option  in  Clemm* s  (1969)  algorithm  for  the  radial  Mathieu  functions, 
and  to  multiply  the  computed  values  by  a  factor  /n/2. 

The  power  series  in  (2.22)  was  explicitly  evaluated,  up  to  a 
truncation  order.  Ihis  truncation  order  was  reached  when  its  addition 
changed  the  summation  value,  in  a  relative  sense,  less  than  a  given  toler¬ 
ance  level.  This  tolerance  level  was  10“4. 

As  a  check  on  the  correctness  of  the  solutions,  qualitative 
comparisons  of  R  and  a  results  with  figures  in  Burke  and  Twersky  (1966) 
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were  performed.  As  a  quantitative  check,  the  expression  (2.26)  for  energy 
conservation  was  evaluated.  lhe  integral  over  outgoing  angles  was 
evaluated  using  a  simple  trapezoidal  numerical  integration.  The  average 
resolution  from  0  to  it  was  about  5°,  with  2°  resolution  in  the  vicinity  of 
the  main  scattering  lobe.  The  error  was  found  to  be  on  the  order  of  1%. 
Better  accuracy  might  be  achieved  by  improving  the  angular  resolution  in 
the  numerical  integration,  or  by  decreasing  the  tolerance  level  in  the 
summation  of  (2.22).  Of  course,  both  approaches  would  lead  to  increased 
computational  expense. 


Section  4 
RESULTS 


This  section  presents  results  for  calculations  of  coherent  power 
reflection  coefficient  R,  and  scattering  cross  section  <j.  For  these  cal¬ 
culations  we  make  certain  assumptions  concerning  ice  keel  shapes  and 
distributions,  we  assume  that  keels  have  a  half-width-to-depth  ratio  of 
1.6,  that  they  have  a  random  orientation,  and  that  they  are  spaced  random¬ 
ly  along  a  track.  Consider  the  assumption  that  keels  have  random  orienta¬ 
tion.  Even  if  the  ice  field  contains  keels  of  only  a  single  size,  the 
projected  keel  width  intercepted  by  any  particular  track  will  take  on  a 
range  of  values.  The  average  projected  keel  width  will  be  equal  to  the 
actual  keel  width  times  a  factor  */2.  If  the  keels  are  spaced  randomly, 
then  we  can  attribute  an  average  separation  distance  to  an  ice  keel  field. 

Figure  4.1  shows  reflection  loss  as  a  function  of  grazing  angle, 
for  three  frequencies;  50,  100,  and  200  Hz.  The  keel  depth  is  4.3  m. 
Therefore  the  keel  width  is  13.76  m,  and  the  average  projected  width  is 
21.61  m.  We  assume  an  average  keel  separation  distance  of  92  m.  The  200 
Hz  curve  increases  most  rapidly  with  grazing  angle,  because  its  small 
wavelength  (7.5  m)  is  best  able  to  resolve  the  keel  shape.  The  curve 
reaches  a  peak  at  30°,  and  then  falls  as  shadowing  decreases  and  more  flat 
surface  becomes  illuminated.  Qualitatively  similarly  shaped  curves  were 
presented  by  Burke  and  Twersky  ( 1 966  ) . 

Figures  4.2  -  4.12  show  scattering  cross  section  a  as  a  function 
of  scattered  angle,  for  different  frequencies,  ranging  from  20  Hz  to  3500 
Hz.  Each  figure  shows  four  incoming  grazing  angles;  5°,  10°,  15°,  and 
20°,  signified  by  small  arrows  on  the  left  side  of  each  plot.  The  reader 
should  note  that  the  radial  coordinate  axes  have  different  scales  on  each 
of  the  plots. 


At  20  Hz  and  50  Hz,  a  single  main  scattering  lobe  appears  in 
each  plot.  Secondary  lobes  begin  to  appear  at  100  Hz  and  above.  At  high 
frequencies,  the  secondary  lobe  structures  become  quite  complicated.  The 
main  scattering  lobes  increase  in  magnitude  with  increasing  frequency  and 
with  increasing  grazing  angle.  They  also  become  more  narrow  with  increas¬ 
ing  frequency. 

The  orientations  of  the  main  lobes  decrease  in  angle  with 
increasing  frequency.  Above  about  250  Hz,  the  orientation  angle  also 
increases  with  incoming  grazing  angle.  At  the  highest  frequencies,  the 
lobes  are  oriented  at  angles  nearly  equal  to  their  respective  incoming 
grazing  angles.  In  other  words,  at  high  frequencies,  most  of  the 
scattered  energy  is  directed  into  the  same  angle  as  that  of  coherent 
specular  reflection. 

We  computed  the  surface  scattering  loss  coefficient  SL,  given 
by 

n 

SL  *  1  -  R  -  csc(u)  /  o($,u)d$  «  csc(u)  /  o(<fr,u)d<fr  .  (4.1) 

0 

Note  the  difference  between  the  ranges  of  integration  in  (4.1)  and  (2.26). 
The  range  from  0  to  $  is  the  angular  range  of  propagation  in  the  Arctic 
Ocean.  Energy  directed  below  *  a  15°,  represented  by  the  last  term  in 
(4.1),  is  eventually  absorbed  by  bottom  loss.  Figures  4.13  and  4.14  show 
R  and  1 -SL  as  a  function  of  frequency,  for  grazing  angles  u  =  5°  and  u  = 
10°,  respectively.  These  figures  show  that  at  frequencies  below  200  Hz, 
the  coherent  power  reflection  loss  is  a  good  estimator  of  scattering  loss. 
However,  above  200  Hz  R  and  1-SL  begin  to  diverge.  The  reason  is  that  a 
significant  fraction  of  the  scattered  energy  is  directed  within  the 
angular  range  of  propagation,  from  0°  to  15°.  Therefore,  it  is  necessary 
to  evaluate  the  full  expression  (4.1)  for  an  accurate  estimate  of  scat¬ 
tering  loss  at  high  frequencies. 
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Figures  4.15  -  4.24  show  the  scattering  cross  section  a,  in  a 
manner  analogous  to  Figs.  4.2  -  4.12,  for  a  different  set  of  keel  param¬ 
eters.  The  frequencies  range  from  20  Hz  to  1600  Hz.  The  parameters  for 
this  case  are:  keel  depth  =*  7.7  m,  average  projected  width  -  38.7  m, 

average  separation  =*  62  m.  This  case  gives  results  that  are  qualitatively 
similar  to  those  of  the  first  case.  Quantitatively,  the  magnitude  of  a  is 
greater  for  this  second  case,  because  of  the  larger  size  and  denser 
concentration  of  keel  cylinders. 

The  results  maybe  summarized  as  follows: 

1)  At  low  frequencies,  below  about  200  Hz,  most  of  the  incoher¬ 
ent,  scattered  energy  is  directed  into  large  angles  with 
respect  to  the  horizontal,  and  therefore  exits  from  the 
sound  propagation  channel. 

2)  With  increasing  frequency,  the  orientation  angle  of  the  main 
lobe  of  the  scattering  cross  section  asymptotically 
approaches  the  incident  grazing  angle.  Above  200  Hz,  a 
significant  fraction  of  the  incoherent,  scattered  energy  is 
directed  into  shallow  propagation  angles,  less  than  15°. 

For  realistic  predictions,  acoustic  propagation  models  must 
incorporate  these  effects.  This  was  done,  for  example,  in  an  energy 
transport  model  described  by  Greene  and  Rubenstein  (1986).  This  transport 
model,  in  fact,  includes  azimuthal  as  well  as  vertical  angle  scattering 
cross  sections.  This  approach  allows  a  more  general  three-dimensional 
prediction  of  propagation. 
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COHERENT  REFLECTION  LO! 


figure  4.1. 


Reflection  Coefficient  R  for  Frequencies  50,  100,  and  200  Hz 
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rigor*  4.2.  Scattering  cro«*  section*  at  20  Hs.  Incoming  grazing  angle* 
are  5*,  10*,  15*,  and  20*.  Keel  parameter*!  Width  -  21.61 
a,  Depth  -  4.3  a.  Average  Separation  -  92  a. 
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Figure  4.4.  Aa  in  Figure  4.2,  but  at  100  Hz , 
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Figure  4.4.  (Continued) 
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Figure  4.6.  (Continued) 
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Figure  4.7.  ha  in  Figure  4.2,  but  et  250  Hz 
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Figure  4.8.  (Continued) 
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Pigure  4.13. 


Reflection  coefficient  R  (solid  carve)  end  1-SL  (syebols) 
for  5*  gracing  angle.  The  quantity  1-SL  (SL  is  scattering 
loss)  represents  the  sen  of  coherent  energy  plus  the  propa¬ 
gating  coaiponent  of  incoherent  scattering  energy. 
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Figure  4.14.  Reflection  coefficient  R  (solid  cur re)  and  1-SI*  (symbols) 
for  10*  grazing  angle.  The  quantity  1-SL  (SL  is  scattering 
loss)  represents  the  sue  of  coherent  energy  plus  the  propa¬ 
gating  component  of  incoherent  scattering  energy. 
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Figure  4.16.  As  in  Figure  4.15,  but  at  50  Hz 
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Figure  4.16.  ( 
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Figure  4. 19.  Aa  in  Figure  4.15,  but  at  200  Hz, 


Figure  4.22.  hm  in  Flgura  4.15,  bat  at  400  Be. 
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Figure  4.23. 


As  in  Figure  4.15,  but  at  800  Hz . 
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Figure  4.24*  (Continued) 
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